Properties of excitations in systems with a spinor Bose-Einstein condensate 
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General theory in case of homogenous Bose-Einstein con- 
densed systems with spinor condensate is presented for the 
correlation functions of density and spin fluctuations and for 
the one-particle propagators as well. The random phase ap- 
proximation is investigated and the damping of the modes 
is given in the intermediate temperature region. It is shown 
that the collective and the one-particle excitation spectra do 
not coincide fully. 

PACS numbers: 03.75.Fi, 67.40. Db, 05.30.Jp 



Since the realization of Bose-Einstein condensate in 
vapours of alkali-metal atoms, there has been a tremen- 
dous development of the field concerning experimental 
and theoretical aspects as well. Most of the works have 
considered the gases in magnetic traps, where the spin 
degrees of freedom are frozen and he particles behave as 
scalar ones. The confinement of the 23 Na gas in an op- 
tical trap Q| has opened a new avenue of investigations 
where the spinor nature of the particle field plays a de- 
cisive role |2J-|lO||. The present paper is devoted also to 
this problem considering a homogenous system. 

Two states of systems with spinor condensates have 
been discovered |§-{||, namely the polar and the ferro- 
magnetic ones and the conditions of their appearance 
have been determined. It has also been shown that the 
collective excitations possess different features in the two 
states influenced also by a possible external magnetic 
field P,H,0- The Bogoliubov approximation, which has 
been applied until now, does not take into account the 
interaction between the particles in the condensate and 
the thermal cloud, leading to undamped modes. 

In this paper we start from the thermal Green's func- 
tion formalism |0-H|. First general relationships are 
derived for the density and spin correlation functions by 
expressing them in terms of their proper parts. Below Tq 
the proper parts can be separated into singular and regu- 
lar contributions such that the diagrams contributing to 
former ones split by cutting a single one-particle propa- 
gator. In lowest order the singular diagrams alone lead 
to the Bogoliubov approximation for the modes. Treat- 
ing the regular part also in leading order damping of the 
modes and other qualitatively new features emerge. The 
one-particle propagators are also investigated and it is 
pointed out that the one-particle spectra coincide with 
a part of the collective ones in the symmetry breaking 
phase. Finally the damping of the modes are determined 
in the intermediate temperature region. 



The model used here is the usual one for particles with 
spin 1. The Hamiltonian is 



H 



E 

k 



(e k - ju)e4(k)a r (k) + 
+ l - J2 4,(k 1 )at(k 2 )F; s ' s 'a s (k 3 )a s Kk4), (1) 

ki+k 2 = 
=k 3 +k 4 

where ek = h 2 k 2 /(2M) (M is the mass of the atom), 
fi stands for the chemical potential and the two particle 
interaction is given by 



V 



C n rs r > s > -f~ C s ^.f j rs yr ) 7 



(2) 



Here the constants c„ and c s are related to the s-wave 
scattering lengths ao and ai in the total spin zero and 
two channel, respectively, by 3c„ = 4irh 2 a n /M and 2c s = 
4irh 2 a s /M, where a n = a + 2a^ and a s = 2(a 2 — ao)/3 
0-0- We assume that the system has unit volume. 
a r (k) and a£(k) are annihilation and creation operators, 
respectively, for particles with momentum k and spin 
component r. The representation r = 1,0,-1 will be 
used, where the spin operators are F z = diag(l,0, — 1), 
{F±) rs = (F x ) rs ± i{Fy) rs = \/2((5 ± i, r <5o, s + b~o, r 5±i,s)- 
In Eq. (fy) and from now on summation over repeated 
indices is understood. 

To take into account the condensation the canonical 
transformation 



Or(k) = & r (k) + 4,OV^oCr 



(3) 



is introduced, where Nq denotes the number of particles 
in the zero momentum state of spin projection r, forming 
the condensate and Q is a normalized spinor. Nq is to be 
determined self-consistently. 

As mentioned before there are two states according to 
the sign of the constant c s : The polar state (c s > 0), 
where (F) = 0, (C P )r — S r .o, and the ferromagnetic state 
(c. <0), where (F) 2 = 1, (( F ) r = 6 rA . 

We will deal with the following correlation functions 



(4) 
(5) 
(6) 
(7) 
(8) 



D nn (k, T ) = (T T [n(k,T)n(-k,0)] 

D zz (k, T ) = (r T [f z {k,T)f z (-k,o)] 

D nz (k )T ) = (T r [n(k,r)J^(-k,0)] 
£>±±(k,r) - (T 7 .[jr ± (k J r)^(-k,0)] 
D QQ (k,r) = (T T [jr Q (k,T)^(-k,0)] 



Here n(k) = J2„al(k + q)a r (q) is the density operator, 
T z (k) = E q 4(k + q)(^) rs « s (q), ^±(k) = E q 4(k + 
q)(F±) rs a s (q), T Q {k) = E q <4( k + q)(^+)r S a s (q) are 
the appropriate spin-density operators. The averaging is 
made over the grand canonical ensemble, r is the imagi- 
nary time and T T is the r ordering operator pj| ] . 

We introduce a propagator from what the above ones 
can be easily evaluated as 

D s r r s ,(k,r) = (T T [<7 rs (k,r)<v r '(-k,0)]Y (9) 



where a rs (k) = E q <4( k + q)a s (q)- 

The propagators are periodic in r with period /3S, in 
their Fourier-series the occurring frequencies are the Bose 
discrete Matsubara frequencies Lo n = ^§§- Q ■ The prop- 
agators DfJ s ,(k, iuj n ) can be written in terms of their 
proper parts n^ s /(k, iuj n ) (which do not split by cutting 
a single interaction line): 



D s r 7 s ,(k,iu n ) = U°7 s ,{k,iu> n ) + 
+W a l(k,iu Jn )V < %DX,(k,iu> n ). 



(10) 



Space reflection and time inversion symmetry yields 
DfJ s , — D r sr s . Spin conservation means that r — s = 
r' — s'. The same symmetry relations hold for the proper 
parts. As a result Eq. dlG) separates into independent 
equations for specified spin transfers. 

The first one (corresponding to zero spin transfer) 
reads as °£ = °g + °g ■ °£ ■ °J2, where the following 
notation is introduced 



>D = 



D ++ 

D 00 

D77 



D ++ 
D 00 



Dn 




(11) 



and equivalently for °IL Here and in the following +, 0, — 
is used instead of +1, 0, — 1 for the spin components. The 
interaction matrix is °C_ = c n ■ £ o £ + c s ■ £ o £ where 

£ = (l,l,l) T and£ = (1,0,— 1) T . One can easily verify 
that after solving the matrix equation for °D_ the corre- 
lation functions defined in (0)-(p[) can be written in the 
following form: D nn =g-°R-£ L , D zz = g ■ °R ■ ^ and 
D nz = £ • °D_ ■ £, ■ Since °_D is a symmetric matrix (due 
to the above discussed symmetry relations) D nz = D zn . 
The second equation (corresponding to plus one spin 
transfer) is + D = +11 + +11 • + C ■ + D_, where 



/ D o+ D o+ 



(12) 



and its proper part matrix + n is similarly defined. The 
interaction here can be written in the form + C_ = c s -x°Xj 
where x = (1j 1) T - Furthermore D ++ = x T ■ + D_ ■ %. 

The equation for the negative spin transfer (for ~ D) 
can be calculated in the same way. The equation cor- 
responding to the plus two spin transfer is trivial, since 



the interaction matrix becomes zero, consequently this 
propagator is already proper: -D++ = n^+. 

The excitation energies and their widths can be found 
by analytic continuation on the omega plane in the 
usual way pd|| and are determined by the condition 

det(J — +n • +C_) = 0, where / is the unit matrix of 
rank three and two in the case of zero and plus one spin 
transfer, respectively. The simplest possible approxima- 
tion is to keep only the leading singular contributions for 
the proper parts. Such an approach is equivalent to the 
Bogoliubov approximation. The next simplest calcula- 
tion is to keep also the leading regular diagrams as well, 
giving the contribution Ho • /. Here n is the contribution 
of the bubble diagram 



n (k,iw„) 



d 3 q n°(k + q)-n°(k) 
(27r) 3 hiuj n - (e k+q - e k ) ' 



(13) 



where n°(k) is the Bose-Einstein distribution function for 
particles of energy e k . 

The calculation in this random phase approximation 
can be carried out separately for the two states. 

In the polar state the leading singular contributions are 



n p o 
o 



h 2^ = 



n? 



(14a) 
(14b) 



In Eq. @ U s = 2N h- 2 e k /[{nu n ) 2 - h- 2 el], U ± = 
N h- 1 /{±iLj n -h- 1 e ]i ). 

Using °g = n • °L+°R B the determinant det(I-°g- 
°£) = (1 - 2c s n )(l - 3c n n - c n n s ). These spectra 
separate into two parts, corresponding to the density and 
the spin density mode described by 



Dnn 



3n + u s 



1 - 3c n n - c n U s 



and 



D zz = 



2U 



i - 2c„n, 



o 



(15) 



(16) 



respectively. 

The density and the spin fluctuations are not coupled 
in this state (D nz = 0) and the spin fluctuations are 
created only by the thermal cloud. 

For the spin wave mode det(I— + n • + C_) = 1 — 2cTIo — 
c s Us and D ++ = \ T • + D_ ■ x results in 



D ++ = 



2n + n s 



1 - 2c s n - c s U s 



(17) 



The same result holds for the D spin-wave mode. 

The quadrupolar spin-density fluctuation propagator is 
simple since Dqq = 4D7_ _ = 4Uo, which means that 
this belongs also only to the non-condensate. Keep- 
ing only the singular contributions II the poles of D nn 



and D ++ are at frequencies iof m — h 1 \/ek(ek + 2N c n ), 
u+ + = fi, _1 -\/ek(ek + 2iVoc s ), respectively, in agreement 
with the modes known from previous calculations |3|j] . 

We do not go into any detail concerning the ferro- 
magnetic state and give the result only for the density 
fluctuations. The determinant A = det(| - °g • °0 = 
(1 - 3c„n - CnIIs)(l - 2c s n - c s U s ) - c„c s n| does not 
factorize into two parts, which means that the density 
and the spin density autocorrelation functions (defined 
in Eqs. (H) and (0)) are coupled, so the cross correlation 
function (defined in Eq. (H)) is not zero either. One can 
find 



Dr. 



D zz = 



3n (i-2c,n ) + n s (i-5c s n ) 

A 

2n (i - 3c„n ) + n s (i - 5c„n ) 



Unz — *-* zn 



ns 

A ' 



(18a) 
(18b) 
(18c) 



where A is the above given determinant. 

We note that the results of the Bogoliubov approx- 
imation can be simply regained here as well by set- 
ting IIo to zero, leading to the eigenfrequency u^ = 
^ _1 \/ e k[ek + 2No(c n + c s )] in agreement with J3), Q. 

Before further discussing the spectra of the correlation 
functions (0)— (H) we present the structure of the thermal 
Green's functions defined by 



^(k,T)=-(T T [^(k,r)6f(k,0)] 



(19) 



The Roman indices are standing for the spin compo- 
nent, while the Greek ones (taken as +1, -1) are for dis- 
tinguishing between the normal and anomalous Green's 
functions fO-H], using the notation bj.(k) = b r (k), 

b-\k) = b]r-kj. 

The generalized Dyson-Beliaev equation for the 



Green's function reads as G™g 



Qrs 



i n.T8 v s s n s s 



where S^ is the self energy, the sum of those diagrams 
with two external points which can not be split by cut- 
ting a single particle line. The free propagator in our 
case is G™ a/3 {k, iu> n ) = 6 rs 5 al3 /(aiu) n - h~ l eu). 

Time reversal and space reflection invariance mean 
that C-(k,iw„) = Gf a {k,iw n ) = G r J a ^ (k,-ico n ). 
Spin preservation requires for Q r \ that ar — f3s — (a — 
f3)a, where a is the spin of a condensate particle. This 
structure can be symbolized for the polar state by 



G X X 



■ 











■ 











■ 






and for the ferromagnetic state by 



r>rs 



■ 
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■ 



r>rs 

y~xx 
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where a black box refers to a nonzero element. The same 
structure holds for the self energies as well. 

The dielectric formalism worked out for scalar Bose 
particles |13|-[l9| can be generalized for the present case. 
In this framework the expressions for the correlation 
functions D nn , -D++ in the polar state and D nn (D zz , 
D nz ) in the ferromagnetic state can be rearranged so 
that they have the same denominator as the propagators 
0?? (£i_i), GH (£1-1) m the polar state and Gil (0^) 
in the ferromagnetic state, respectively. (We listed only 
those relevant to our considerations here.) Consequently 
the above pairs share the same excitation spectra. 

We are not going to present the general theory here, 
instead we demonstrate this property in the RPA. The 
interaction potential defined in Eq. (|2|) is renormalized 
to 



W£'(k,iw„) 



1 - 3cnILo(k,iu>n 



-0 rs r > s ' \ 



1 - 2c s U (k,iuj n ) 



(F) rs (F) r 



(20) 



The diagrams contributing to the self-energies are drawn 
in Fig 0. 

For the polar state the resulting Green's functions are 



(i - cii){iuj n + Tr^fk) + fi-ijVoc 



C++ 

y Xl - i, , . ,2 t-2,,2 



[{iuj n ) 2 - h- 2 el]{\ - en) - 2h- 1 e k N c 

where c = c n , II = 3Ho and c = c s , II = 2Hq for Gii and 
C++, respectively. 

For the discussed mode in the ferromagnetic state the 
Green's function is 



C++ - 



h(ihw n + e k )(l - 3c T JI )(l - 2c s IIo) + HN p 
[(ihiu n ) 2 + e 2](l - 3c.no) (1 - 2c s n ) + 2N e kP ' 

(22) 

+ c s — 5c n c s Tlo is introduced 



Here the abbreviation p = 
to simplify the equation. 

Comparison of Eq. (|2l]) with ( |l5| ) and ([17]) and further- 
more (|23) with (U^) illustrate the statement above. Note 
that D^z and Dqq do not have counterparts among the 
Green's functions, which means that the collective and 
one-particle spectra do not coincide completely. 

Now we turn our attention evaluating the spectra. It 
is useful to introduce characteristic lengths, namely the 
thermal wavelength A = ^/(2A/fcsT) 1 / 2 , the correla- 
tion length £' = MksT/(4:Trh 2 No) governing the critical 
fluctuations and the mean field (Bogoliubov) coherence 
lengths ^ s = h/iiMNoc^s) 1 / 2 . 

We consider the intermediate temperature region, 
where £ B ^> A, £' and the contribution of the bubble 
graph can be treated as a perturbation. To leading or- 
der shift in the energy can be neglected and the main 
effect is the damping of the eigenmodes due to their in- 
teraction with the thermal cloud represented by the bub- 
ble graph. This is a Landau-type damping which is the 
dominating process at such temperatures, no has been 



analyzed in detail in J14j . In the temperature region con- 
sidered and for linear dispersion Ilo can be approximated 
by i[M / {Ai:h 2 k\ 2 )] ■ h\[(hw - e k )/(fuj + e k )]. The spectra 
can be written in the form uj—ij, where w is the frequency 
in the Bogoliubov approximation for small k limit and 7 
represents the damping if u) is not zero. In the polar state 
hL/"L) = (3f«) and ( 7 £+K + ) = W/%f)- 
Furthermore D zz has a pole at ujf z — 0, j zz = 27^ + . In 
the ferromagnetic state 7^ = j£ z = 7^ and 7^„/w^ n = 



In the experiments made with sodium atoms [l), the 
scattering lengths are ao — (46 ± 5)as, 0,2 = (52 ± 5)as 
(see e.g. |3J), where as is the Bohr radius, which means 
that c n > i.e. the polar state is realized in equilibrium 
and in the absence of an external magnetic field. For 
an estimate the temperature can be taken to be of the 
order of 10~ 7 K and the condensate density is about of 
3-10~ 20 m~ 3 . With these parameters 7^/ ui^ w 1,2-10 -2 
and 7^/w^ « 1,6- 10~ 3 , showing that the spin wave is 
much less damped than the density fluctuations. 

Finally some remarks are appropriate about the rele- 
vance of calculations for homogenous system to the Bose 
gas in a trap, beyond the obvious usefulness namely that 
they can serve as a guide. As far as the damping of the 
modes is concerned one can refer to the fact that the 
simple expression obtained in RPA []l2| , [L4[ has proved to 
be quantitatively a good approximation |l9-G2j in case 
of gases in magnetic traps, even when compared to the 
experimental findings |23|. One may expect a similar 
situation here. Moreover, the results obtained for the 
homogenous system are relevant in the local density ap- 
proximation for the gas in a trap. Namely by exciting 
density perturbations much smaller than the size of the 
condensate the local speed of the quasiparticle (zeroth 
sound) has been measured and its density dependence 
has been found in agreement with the bulk result p4| . 
One can hope that such experiments can be extended to 
spin waves. 
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FIG. 1. Self energy contributions in RPA, the solid line 
denotes a one-particle free propagator, a double-dashed-line 
denotes the renormalized interaction and the circle is repre- 
senting a conensate particle. 
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